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Abstract 

The mass and the total decay width of the W boson are measured with the L3 
detector at the LEP e^e~ colhder using W-boson pairs produced in 0.7 fb~^ of data 
collected at centre-of-mass energies between 161 and 209 GeV. Combining semi- 
leptonic and fully-hadronic final states, the mass and the width of the W boson are 
determined to be 

mw = 80.270 ± 0.046 ± 0.031 GeV and 
Fw = 2.18 ±0.11 ±0.09 GeV , 

where the first uncertainty is statistical and the second systematic. 
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1 Introduction 



The mass, mw, and the total decay width, Fw, are fundamental properties of the W boson. 
Their measurement, initially performed at the SppS hadron collider [1], provides important 
information about the Standard Model of electroweak interactions [2]. Together with other 
electroweak parameters such as the Z-boson mass, the effective weak mixing angle and the 
measurement of the top-quark mass [3], the precise determination of m-^ allows a thorough 
test of the Standard Model at the quantum loop level as well as constraining the mass of the 
Higgs boson [4]. 

In e"^e~ collisions, W bosons are produced singly or in pairs. At centre-of-mass energies, y/s, 
exceeding 2m-w, W-boson pair production, e+e~ W"'"W~, dominates. The pair-production 
cross section at threshold is sensitive to m-w- Therefore, at LEP m-yv was first derived from 
cross section measurements [5,6]. At centre-of-mass energies well above production threshold, 
W bosons are directly reconstructed and the effective mass of the decay products is used to 
determine mw [7,8]. The mass distribution of the W bosons is analysed and mw and Fw are 
determined by comparing samples of Monte Carlo events to data. A reweighting procedure is 
applied to obtain Monte Carlo samples corresponding to different values of mw and Fw- 

In the following, mw and Fw are defined such that the denominator of the W-boson 
propagator, (m^ — mw) + im^Fw/^w, models the mass-dependent width of the W boson. 
The analysis presented here is based on a data sample collected with the L3 detector [9] at 
y/s — 189 — 209 CeV, corresponding to a total integrated luminosity of 629 pb~^. These results 
are combined with previous L3 measurements at lower centre-of-mass energies [5, 7] yielding 
final results on mw and Fw based on the complete data sample of 706 pb~^ collected by the L3 
experiment at i/s = 161 - 209 GeV. Other measurements of mw at LEP and the TEVATRON 
are described in References 8 and 10, respectively. The indirect determination of mw from 
electroweak precision data is presented in Reference 4. 

2 Data sample 

W bosons decay into hadrons, mostly through ud or cs, or leptons, W~ — > i^i^e, where 

i denotes an electron, muon or tau lepton. Charge-conjugate states are understood to be 
included throughout this article. In the following, these final states are denoted as qq and iu, 
or, in general, //, for both W"*" and W~ decays. W-boson pair production yields three classes 
of events: the fully-leptonic, iuiu, the semi-leptonic, qqiu, and the fully-hadronic, qqqq, final 
states. Due to the presence of more than one neutrino in the iuiu final state, the effective masses 
of the W bosons cannot be directly reconstructed from their decay products and this decay 
channel is not further considered here. Visible final-state fermions are reconstructed in each 
event. Electrons and muons from W-boson or r-lepton decays are measured in the calorimeters 
and in the tracking system. Hadronically-decaying r-leptons are identified as narrow, low- 
multiplicity, jets. Jets originating from quarks are reconstructed by combining information 
from calorimetric clusters and associated tracks into jets using the DURHAM algorithm [11]. 

The data analysed correspond to seven average values of ^/s, listed in Table 1. The selection 
of W-boson pair- production events is described in Reference 12. The selection of the qqeu and 
qqjUi^ final states requires an identified high-energy electron or muon, respectively. The qqru 
final state is characterised by a low-energy isolated electron or muon or by the reconstruction 
of a narrow jet. For all semi-leptonic final states, missing momentum due to the neutrino is 
required and the jet-jet mass has to be compatible with mw- The selection of the qqqq final 
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state requires events with high multiphcity, small missing energy and a four-jet topology. To 
reject quark-pair production with additional jets originating from radiated gluons an artificial 
neural network is trained using discriminating variables such as the jet energies, broadenings 
and angles, the event spherocity, the jet multiplicity and the DURHAM jet-resolution parameter 
for which the event topology changes from three to four jets, 7/34. Only events with high neural- 
network output are retained for further analysis. The numbers of selected W-boson pairs are 
detailed in Table 1. 

3 Monte Carlo simulation 

The KANDY [13] Monte Carlo generator is used to model four-fermion production, including 
both W-boson production and background processes. The RACOONWW [14] program is used 
as a cross check and to estimate possible systematic uncertainties due to the modelling of photon 
radiation. Additional background contribution from fermion-pair production, dominated by the 
e+e" qq process, is simulated using the KK2F [15] event generator. Monte Carlo events are 
generated at the seven average ^/s values listed in Table 1. Effects from the spread of centre- 
of-mass energies within the individual energy points are found to be negligible. The expected 
number of background events is listed in Table 1. 

The hadronisation process is modelled with the PYTHIA [16] program, while the HER- 
WIG [17] and ARIADNE [18] programs are used to assess systematic uncertainties. These 
Monte Carlo programs are tuned to describe hadronic Z-boson decays recorded at the Z reso- 
nance [19]. In the case of W-boson pair production, a dedicated parameter set, tuned only on 
Z-boson decays into light-quarks (u,d,c,s) is used. 

Bose-Einstein correlations (BEC) [20] in W-boson decays are simulated using the BE32 
model [21] implemented in PYTHIA. Only BEC between hadrons originating from the same 
W boson are taken into account, as suggested by our measurements [22]. Colour- Reconnection 
(CR) effects [23] in the qqqq final state would alter the colour fiow between the W bosons. 
In accordance with our measurement [24], these are not implemented in the Monte Carlo 
simulation. However, both the effect of BEC between hadrons originating from different W 
decays and that of CR between W bosons are considered as possible systematic uncertainties. 

The response of the L3 detector is modelled with the GEANT [25] program which includes 
the effects of energy loss, multiple scattering and showering in the detector material. Hadronic 
showers are simulated with the GHEISHA [26] program. Time-dependent detector efficiencies, 
as monitored during data taking, are included in the simulation. 

4 Event reconstruction 

In the qqez/, qq/iz/ and qqqq channels a kinematic fit is applied to improve the resolution of 
the measured energies, Ef, momenta, p/, polar, 9f, and azimuthal, 0/, angles of the visible 
fermions. Four-momentum conservation and other constraints, as detailed below, are imposed. 
The measured quantities are varied within their resolution to satisfy these constraints. The 
resolution of each individually-measured object depends on details of the reconstruction, such 
as the detector region or the energy scale. The average resolutions of Ef, 6f, and 0j for 
electrons, muons and hadronic jets, as determined by Monte Carlo simulation, are given in 
Table 2. These values agree with the resolutions derived from calibration data collected at the 
Z resonance within the statistical accuracy of the test. 
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In all events, four-momentum conservation is required which, in the case of the qqev and 
qqfiu final states, determines the momentum and the direction of the neutrino. For hadronic 
jets, the velocity (3f = Pf/Ef is fixed to its measured value, as many systematic effects cancel in 
this ratio. When imposing energy conservation, the ^/s value determined for each event by the 
LEP Energy Working Group [27] is used. Events collected during the manipulation of the LEP 
beams, for which no precise calibration of the LEP energy is available, are excluded from the 
analysis. Energy conservation results in a one-constraint (IC) kinematic fit for qqez/ and qqfiu 
final states and a four- constraint (4C) kinematic fit for the qqqq channel. In general, fermion 
angles are better measured than energies and momenta. Therefore, the kinematic fit improves 
more the determination of the latter. The improvement in the resolution of the average value 
of the two reconstructed W-boson masses is shown in Table 3. 

In the IC and 4C fits, the masses of the two W bosons are determined separately. The 
mass resolution is further improved by the additional constraint of requiring these masses to be 
equal within the width of the W boson, fixed as 2.1 GeV. This numerical value does not bias 
the resulting fits. This procedure results in a two-constraint (2C) fit of qqeu and qqjiiy events 
and a five-constraint (5C) fit of qqqq events. For qqeu and qq/iv events, both the mass of the 
hadronically-decaying W boson obtained in the IC fit, rriic, and the average mass obtained in 
the 2C fit, m2c, are used in the mass extraction, which is described in the following section. 
Similarly, in the qqqq channel, the average masses of the 4C fit, m^c, and the 5C fit, m^c, are 
used. 

The qqru final state contains at least two neutrinos and only the W boson decaying into 
hadrons is used in the mass reconstruction. The energies of the two hadronic jets are rescaled 
by a common factor such that the sum of their energies equals \/s/2, effectively imposing an 
equal-mass constraint on the two W bosons. Use of the mass of the hadronic system after this 
rescaling, mresc; improves the resolution of the W-boson mass reconstruction by more than a 
factor of two. 

The improvement of the mass resolution due to the kinematic fit is shown in Figure 1. 
The average mass resolutions before and after the kinematic fits or the energy rescaling is 
summarised in Table 3 for all final states. Only the better-measured quantities m2c and m^c 
are used to determine Fw 

W-boson pair production is frequently accompanied by photon radiation. Photons near to 
a final-state fermion are mainly due to final-state radiation (FSR). In qqez/ events, photons 
close to the electron are automatically included into the measurement of the electromagnetic 
cluster. In qqjiu events, the cluster closest to the muon direction is assumed to originate from 
the ionisation energy loss of the muon in the calorimeters and is taken out of the event. In qqru 
events, the photon clusters are combined in the tau jet by the jet-reconstruction cone algorithm. 
Hard photons with energies greater than 5 GeV and outside a cone of 5° half-opening angle 
around the lepton are detected in 5% of the qqez/ events and in 2.5% of the qquz/ events. They 
are taken into account by the kinematic fit, but not incorporated in the mass reconstruction, 
as they are mainly due to initial-state radiation (ISR). In all other cases, the detected photons 
are assigned to the jets during the clustering process. For photons emitted along the beam 
direction, and therefore undetected, the analysis relies on the Monte Carlo simulation. 

Systematic uncertainties arise in the qqqq channel due to potential effects of CR between 
the jets from different W bosons. To reduce these effects, clusters with an energy below 2 GeV 
are removed from the original jets obtained by the jet clustering process, as discussed in Sec- 
tion 6.9. The jet energies and momenta are re-scaled with an equal scale factor in order to 
obtain the original jet energy. Only the jet directions and the jet masses are affected by this 
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procedure and the energy resolution of the jets is preserved as illustrated in Table 2. On the 
other hand, the angular resolution of the jets is worsened leading to a degradation of the W- 
boson mass resolution by about 20%. The resulting increase of the statistical uncertainty on 
mw is overcompensated by a reduction of the systematic uncertainty, leading to a lower total 
uncertainty, as discussed in Section 6.9. 

The mass resolution in the fuUy-hadronic final state is improved by taking into account gluon 
radiation from quarks. The DURHAM jct-rcsolution parameter for which the event topology 
changes from four to five jets, 1/45, is used to separate events with and without gluon radiation. 
Those with log 1/45 > —6.2 are treated as four-jet events and the remaining as five-jet events. 

The four or five jets must be associated with the two W bosons. In the case of four jets, 
all three combinations are considered. Five jets can be paired in ten different ways. Monte 
Carlo studies show that in five-jet events only some combinations have a high probability to 
be correct. These are the ones in which the W boson that decays without hard-gluon radiation 
is formed by the highest-energy jet and any other jet or by the second highest-energy jet and 
any other jet except the lowest-energy jet. Only these six combinations are considered. The 
three pairings with the highest kinematic-fit probability are retained. They are ordered by 
their fit probability and treated as separate samples. Pairings where the fit did not converge 
are rejected. This criterion removes 5% of the events. 

The four-jet and five-jet samples, of about equal size, are treated in separate mass fits since 
their mass resolutions are different by about 30% as shown in Table 3. Due to the overall 
improvement in mass resolution, the statistical uncertainty of mw, as determined in the fuUy- 
hadronic channel, is reduced by 6%. 

The mass spectra after the kinematic fit for the better-measured mass variable mi are 
shown in Figure 2 for the semi-leptonic final states and the best pairing in the fuUy-hadronic 
final state. Figure 3 presents the sum of the semi-leptonic distributions, while Figure 4 shows 
the sum of all four distributions. 



5 Extraction of mw and Fw 

A maximum-likelihood method is used to extract m-w and Fw from the reconstructed masses 
of each event. The extraction of mw and Fw is done separately for each of the four final states, 
qqeu, qqjiiy, qqru and qqqq, and the seven average values of ^/s. For each of these 28 event 
samples a hkelihood function, L(m^,r^), is constructed from the product of the individual 
likelihoods. These are evaluated for each mass reconstruction, i, performed for a given semi- 
leptonic event or a given pairing of the four- and five-jet samples of the fully-hadronic final state. 
Correlations between the reconstructed masses from different pairings are found to be negligible. 
The individual likelihoods are calculated from the normalised differential cross sections in terms 
of the reconstructed masses, mi and m2, 

^' ^ y dmi dm2 \dmi dm2 J 

/(mfit,Ffit)(j(mfit,Ffit) + aback 

where a and (Jback are the accepted signal and background cross sections of the corresponding 
final state. As summarised in Table 3, the masses are chosen as mi = m2c and m2 = mic for 
qq^u and qqeu final states and mi = msc and m2 = m4c for fully-hadronic events. For qqrz/ 
events the doubly-differential cross section is reduced to a singly-differential one and only the 
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rescaled mass of the hadronic system is used, mi = mresc- The normahsation factor f{m^, F^) 
is calculated such that the sum of the accepted background and the reweighted signal cross 
section reproduces the measured cross section. This procedure determines mw and Fw solely 
from the shapes of the mass distributions. In the fits to determine mw, the Standard Model 
relation Fw = 3Gf^w(1 + 2Q;s/37r)/(2v^7r) is imposed [28]. When Fw is extracted, mw and 
Fw are treated as independent quantities and the doubly-differential cross section is reduced to 
a single one, since only the better- measured quantity mi is used for the determination of Fw 
The total and differential cross sections of signal and background accepted by the event 
selection are determined using Monte Carlo simulations. Except for single-W production, the 
background cross sections are independent of mw and Fw- The signal Monte Carlo simulation, 
which is originally generated using a particular value of the W-boson mass, m^, and width, 
is modified in a reweighting procedure to represent a different W-boson mass, m^, and 
width, F|^. Each signal Monte Carlo event, j, is given a new weight, Rj, defined by the ratio 

R rm«* F^* m- F-) - '^1' ^w, r^)]' 

where A4 is the matrix element of the four-fermion final state under consideration. The matrix 
elements are calculated for the generated four-vectors of the four fermions, p""^ "'^ using the 
program EXCALIBUR [29]. Since this program is based on four-fermion final states without 
additional photons, the momentum sum of any ISR photons present in the Monte Carlo events, 
k'j , is taken into account by boosting the four fermions into the rest frame of the event after 
the ISR photon emission. Photons not emitted in the initial state are recombined with the 
closest final-state fermion. It was verified that this method is equivalent to using the KANDY 
program, which simulates photon radiation in the event generation. As a cross check of the 
matrix element reweighting the event weights are evaluated from a Breit-Wigner function. 
Consistent results are observed. 

The total accepted signal cross section for a given set of parameters, and T^, is 

i 

where o"^*^" denotes the cross section corresponding to the total Monte Carlo sample containing 
^gen gygnts. The sum extends over all Monte Carlo events, j, accepted by the event selection. 
The total background cross section is 

^gen 

CThack = X] ATgen ' -^b^k,; > (4) 
; -''back,/ 

where, for each background process I with generated cross section cr^ack,/) ^h^k,i -^back,/ 
the numbers of generated and accepted Monte Carlo events, respectively. 

To determine the accepted differential cross section for a given data event, i, the box 
method [30] is applied. When combined with the reweighting procedure, this method takes 
into account detector and selection effects, efficiencies and purities which depend on mw and 
Fw and correlations between the input masses mi and m2. The accepted differential cross 
section is determined by averaging signal Monte Carlo events inside a two-dimensional mass 
domain, Qi, centred around each data event. To take the different resolutions of mi and m2 
into account, these masses are rescaled by their resolutions cxi and a2 whose averages are given 
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in Table 3. The size of each domain is hmited by requiring a sufficient number of Monte Carlo 

events in the domain. In the rescaled parameter space, the distance, dij, of each Monte Carlo 
event, j, from the given data event, i, with reconstructed masses (mi)^ and (777.2)2 is calculated 
from 



dij 



\ 




(5) 



(71 

and the 400 closest Monte Carlo events are retained. 

The most distant Monte Carlo event, jmax, determines the mass intervals around the data 
event, (Smi)i = |(777i)j-(r77i)j^^^| and {Sm2)i = 1(7772)^- (7772)^^^1, which vary between 200 MeV 
and 600 MeV. After summing the weights Rj of all Monte Carlo events associated to the mass 
domain Qi around the considered data event, the differential cross section of the signal processes 
is given by 



^ d777id7772 Ji TT{5mi)i{5m2)i N^"'' 



i^,(777«^,^^,7n|^^^f^^). (6) 



For the background Monte Carlo simulation, the same domain size as for the signal is chosen 
and the differential distribution of the background is determined from the number of selected 
background Monte Carlo events, {N^^^)i, associated with a given data event: 

/ ]2 \ 1 gen 

/ CI 0"back \ _ J- CTback / Arsel \ /7\ 



^dmi (11712) A 7r{5mi)i{5m2)i N^l^-k 



Onc-dimcnsional boxes in the 77ii space arc constructed for the determination of Fw The 
size of each bin is defined by requiring at least 200, but not more than 1000, Monte Carlo 
events. The bin size is at most ±250 MeV around (7771 )j and decreases to about ±30 MeV 
around the peak of the mass spectrum. For the background Monte Carlo simulation, the bin 
size is chosen as ±1 GeV around (777i)i. 



6 Systematic uncertainties 

The systematic uncertainties on 777w and Fw are summarised in Tables 4 and 5, respectively. 
They arise from various sources correlated or un-correlated between the final states and between 
the various i/i values. The different sources of the systematic uncertainty are detailed in the 
following subsections and their correlations arc discussed in Section 7. 

Systematic uncertainties are assessed by determining At/Iw and AF-yv which are defined 
as the changes of the 777w and Fw results if alternative detector calibrations, Monte Carlo 
simulations or reconstruction procedures are used. Two methods are used for the evaluation of 
Ar77w In the cases where the effect of an alternative Monte Carlo simulation is studied, the 
usual mass fit is used, but the data events are replaced by a high-statistics sample from the 
alternative simulation. The fit result, m^, is compared to the nominal W-boson mass common 
to both Monte Carlo samples, 777|^'\ deriving Am-w = — m^. A similar procedure is used 
to derive AFw- In the cases where the agreement between data and simulation is analysed, the 
shift of the reconstructed mass is calculated for each data and Monte Carlo event. The average 
mass shifts of the data and Monte Carlo distributions are compared to determine A777w 
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6.1 Calibration of -y/s 

The value of ^/s is used as a constraint in the kinematic fit. A variation of ^ would imply 
a shift of the reconstructed masses. The relative uncertainty on mw is the same as that on 
y's, while Fw is less affected. This is verified by comparing simulated event samples in which 
the value used in the kinematic fit is systematically changed. The dependences of mw 
and Fw on ^/s are taken into account using the LEP energy determined for the exact time 
each W-boson pair was recorded. The LEP beam energy is known with an accuracy between 
10 and 20 MeV [27]. The complete error matrix from Reference 27 is used to determine the 
uncertainties on mw and Fw given in Tables 4 and 5, which are correlated between all final 
states. 

As a cross check of the \fs calibration, events from the e"'"e~ Z7 process with hard ISR 
were used to measure the mass of the Z boson [31]. The Z-boson mass, mz, was determined to 
be 91.272 ± 0.032 ± 0.033 GeV, where the first uncertainty is statistical and the second system- 
atic, in agreement with the value measured at the Z resonance [32], mi = 91.190 ± 0.003 GeV. 
Assuming this value of mz, the method determines the average to be 175 ± 68 ± 68 MeV 
lower than the value given by the LEP energy calibration, but consistent within the experi- 
mental uncertainty. 

The intrinsic energy spread of the beams causes a ^/s distribution of the individual events 
with a Gaussian width of 240 MeV. To assess this effect, the ^/s constraint in the Monte Carlo 
events is varied by the same amount. The changes of mw and Fw are negligible. 

6.2 Lepton measurement 

The measurement of the lepton energy in qqcz/ and qqjj.u events affects the mass reconstruction, 
while in the qqrz/ final state it is solely based on the measurement of the jets. Control samples 
of events from the e^e~ — > £~^£~ process are selected in calibration runs at the Z resonance 
and are used to cross check lepton reconstruction. The absolute energy scales for electrons and 
muons are known with a precision of 50 MeV. Varying the lepton energy scale by this amount 
and increasing the lepton energy resolution in the simulation by 25% of the value measured 
with Z-resonance data, results in the changes of mw and Fw detailed in Table 6. Effects due 
to the determination of the lepton angles are negligible. 

The distributions of the energy of calorimetric clusters around the charged lepton are shown 
in Figure 5a and 5b. These clusters are normally joined to one of the jets. If they are not 
correctly described by the Monte Carlo simulation, this might result in a bias on the value of 
mw- To assess this effect, all clusters within a cone of 5° half-opening angle around the lepton 
are excluded from the jets. No significant effect on mw is observed. 

6.3 Jet measurement 

The measurements of jet energies and directions affect the mass spectra and are a potential 
source of systematic uncertainties on mw and Fw- These uncertainties are assigned by varying 
the jet-energy scale by 50 MeV, smearing the jet energies by 1% and smearing the jet directions 
by 0.5°. The sizes of these variations correspond to the uncertainties estimated from e'^e" — > qq 
events collected in calibration runs at the Z resonance. These variations are applied to the 
Monte Carlo sample taken as reference to extract mw and Fw from the data. The effects on 
mw and Fw are given in Table 6. As expected, the largest effect appears in the qqru channel, 
where only the rescaled jets are used and no additional constraint is applied. 
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Generally, the event primary-vertex is shifted with respect to the geometrical centre of the 
detector. If this shift was left uncorrected, it would imply a systematic distortion of the jet 
angles. The actual position of the primary vertex is measured using data and corrected for in 
the reconstruction procedure. The shift is found to be less than 4 mm along the beam axis and 
0.5 mm in the transverse plane, with an uncertainty of less than 5%. Figure 6a shows for each 
data event the shift of the reconstructed mass due to the vertex correction. Assuming that the 
vertex is displaced within the uncertainty of its determination results in a change of mw of less 
than 1 MeV. 

Deviations of the calorimeter positions from their nominal locations would also cause angular 
distortions. To check the angular measurement of the calorimetric clusters the measurement 
of mw is repeated using only clusters associated with tracks. For each event, the mass is first 
reconstructed using the angular information of the clusters and then from the angles of the 
associated tracks. These measurements are independent. The resulting mass-shift distribution 
is shown in Figure 6b. Combining all final states we obtain a change of mw of — 1 ± 9 MeV 
between both methods, consistent with zero. 

Angle-dependent effects in the energy scale of the calorimeters could lead to an additional 
bias in the measurement of the jet angles. For instance, if forward clusters had a relative bias 
towards lower energy than clusters in the central part of the detector, the direction of the jet 
would be shifted towards the central detector region. This effect is expected to be most evident 
in the qqqq events, which are strongly constrained by the kinematic fit. To assess this effect the 
raw jet energies are compared to the jet energies after the kinematic fit for various polar-angle 
regions of the detector. No significant change in mw is observed if the cluster-energy scale in 
the simulation is changed for each polar-angle region to agree exactly with the data. 

The energy spectrum of the clusters and the energy flow with respect to the jet axis are also 
investigated. They are shown in Figures 7a and 8a, respectively. Figures 7b and 7c present the 
effect on mw when clusters below a given energy cut are removed. The changes of mw stay 
within the statistical uncertainty of the test when varying the energy cut from the default values 
of 100 MeV and 2 GeV for qqiu and qqqq events, respectively. No signiflcant change of mw is 
observed if clusters outside a cone around the jet axis are removed, as shown in Figures 8b and 
8c for cones of half-opening angles from 30° to 180°. 

6.4 Fit procedure 

The flt procedure determines mw and Fw without any bias as long as the Monte Carlo sim- 
ulation correctly describes effects such as photon radiation and detector resolution. This fit 
procedure is tested to high accuracy by fitting large Monte Carlo samples, typically a hundred 
times the size of the data sample. The fits reproduce well the generated values mw^ and 
within the statistical accuracy of the test, over a range of ±500 MeV in mw and ±600 MeV in 

Tw- 
in the fit of mw the number of events per box is varied between 350 and 450, while in the 

fit of Fw the minimal number of events is varied between 150 and 250. In addition, the fit is 

restricted to masses in the range between 70 GeV and 90 GeV. No statistically significant effect 

on mw or Fw is observed for any of these variations. 

The reliability of the uncertainties given by the fit is tested by fitting for each final state 

several hundred small Monte Carlo samples, each the size of the data sample. The width of 

the distribution of the fitted central values agrees well with the mean of the distribution of the 

fit uncertainties. 
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6.5 Background 

Background which is not correctly described by the Monte Carlo simulation, either in the total 
number of events or in their mass distribution, could cause a shift of mw and Fw- For both 
the semi-leptonic and the fuUy-hadronic selections, the four-fermion background is scaled by 
±5%. Additionally, background from the e"'"e~ qq process is scaled by ±5% and the slope 
of its mass spectrum is varied by ±10% over the mass range between 65 GeV and 95 GeV. 

The dominant background in the fuUy-hadronic selection is due to e"'"e~ — > qq events with 
multiple gluon radiation. To better reproduce the four-jet rate observed in hadronic Z decays, 
a reweighting of the e+e" qq Monte Carlo events according to the value of 1/34 is applied 
in our standard mass-extraction procedure [12]. Removing this reweighting changes the total 
background contribution by 12% and shifts mw and Fw by 10 MeV and 80 MeV, respectively. 
Half of the shift is taken as systematic uncertainty. 

The effects on mw and Fw due to the variation of height and slope of the background mass 
spectrum and the uncertainty due to the reweighting of the 1/34 spectrum are summarised in 
Table 7. The individual sources are added in quadrature to yield the systematic uncertainty 
due to the background simulation. 

6.6 Monte Carlo statistics 

The systematic uncertainty due to the limited size of the signal Monte Carlo sample used for the 
box fit is estimated by dividing it into several sub-samples of equal size and using each of them 
to fit the data. The systematic uncertainty due to Monte Carlo statistics is then determined by 
extrapolating the spread of the fit results to the full Monte Carlo sample. The total systematic 
uncertainties on mw and Fw due to limited Monte Carlo statistics are given for each final state 
in Tables 4 and 5, respectively. 

6.7 Photon radiation 

Four-fermion production, including its radiative corrections, is modelled by the KANDY and 
RACOONWW Monte Carlo generators. Both programs use pole expansions [33] for the cal- 
culation of 0{a) corrections. KANDY models ISR using the Yennie-Frautschi-Suura (YFS) 
exponentiation scheme [34] , while FSR is simulated by the program PHOTOS [35] in the case 
of charged leptons and by PYTHIA in the case of quarks. Interference between ISR and 
FSR is neglected. RACOONWW implements the full 0{a) matrix element for the radiative 
four-fermion production, e"'"e~ f f f fl- Higher-order corrections coming from multiple ISR 
photons are implemented using a structure-function ansatz. As the calculations implemented in 
RACOONWW are based on massless fermions, the FSR simulation exhibits a minimum cut-off 
on the photon-fermion angle. 

The radiation of hard and isolated photons is better simulated by RACOONWW which 
implements the complete matrix element of the ffffj final state. On the other hand, soft and 
coUinear photons are not generated, which makes the KANDY approach more appropriate for 
comparison with data. KANDY uses a W propagator with a mass-dependent term containing 
the W width, whereas RACOONWW uses a constant term. Because the definitions differ by 
27 MeV in the position of the W peak [36], the mf^"" input to RACOONWW is chosen 27 MeV 
lower than for KANDY in order to give an identical W-boson lineshape. 

A total of 300 000 Monte Carlo events of the qqeu, qq/iu and qqqq final states are generated 
with the RACOONWW program at ^/s — 207 GeV, including full detector simulation. Events 
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with hard- photon radiation are selected at the generator level using the CAL05 algorithm [37] , 
which recombines soft and collinear photons with the nearest fermion. These events, after 
detector simulation, are used instead of data in the mass fit which relies on KANDY as the 
reference Monte Carlo. The change in mw from the comparison of the programs is derived 
and scaled by the fraction of events with hard-photon radiation, which is of the order of 10%. 
The same effect as observed in the qqjiiy channel is assumed for the qqri/ channel where no 
events were generated. In an additional test, the KANDY events are reweightcd such that they 
represent the 0(q;^) ISR corrections instead of the 0{a'^) calculation. The changes of m-yv and 
Fw resulting from these tests are detailed in Table 8. For each final state they are added in 
quadrature to estimate the systematic uncertainties on mw and Fw due to the modelling of 
photon radiation, given in Tables 4 and 5. 

6.8 Hadronisation 

The hadronisation process is modelled by three different schemes as implemented in the Monte 
Carlo programs PYTHIA, HERWIG and ARIADNE. For the perturbative phase PYTHIA and 
HERWIG simulate a parton shower, while a dipole cascade is produced in ARIADNE. The 
Lund string-hadronisation model is used by PYTHIA and ARIADNE, while HERWIG employs 
a cluster model. A comparison of the mass distributions of the three different models with data 
is shown in Figure 9. Within the statistical accuracy, all three Monte Carlo distributions are 
compatible with the data. 

The results for and Fw presented in this paper are based on the PYTHIA model. 
Systematic effects due to modelling of the hadronisation process are determined by comparison 
with the other two programs. In the mass-extraction fit the data events are replaced by high- 
statistics samples of Monte Carlo events generated with PYTHIA, HERWIG and ARIADNE. 
These Monte Carlo samples consist of events which are identical at the four-fermion level and 
thus differ only in their hadronisation. The changes of mw and Fw due to the use of HERWIG 
or ARIADNE are listed in Table 9. For mw, HERWIG and ARIADNE reproduce the PYTHIA 
results within the statistical uncertainty, except in the qqru channel. This is mainly caused 
by the misassignment of energy deposits from the remainder of the event to the tau-lepton jet 
which is based on a cone definition. This effect, altering the reconstruction of the jets and 
therefore mw, strongly depends on the choice of the hadronisation model. For Fw, ARIADNE 
is in good agreement with PYTHIA, while HERWIG shows significant differences, especially 
for semi-leptonic final states. 

The four-momenta of calorimetric clusters which are used to form hadronic jets are calcu- 
lated using the energy and angle measurements and assuming their masses to be either zero 
or the pion mass. However, kaons and protons are frequently produced resulting in a shift of 
the jet masses. This shift is automatically corrected in the mass-extraction fit which uses the 
Monte Carlo simulation containing these hadrons. If the simulation predicts different multiplic- 
ities for these heavier hadrons than is present in data, systematic effects on the measurement 
of mw and Fw are expected. In order to assess this systematic effect, the mean number of 
charged kaons and protons produced in the W-boson decays of our simulation is compared to 
measurements [38] and found to be in agreement. The shifts Amw and AFw are calculated 
with Monte Carlo events reweighted such that the mean kaon and proton multiplicities agree 
exactly with the measured values. It is checked that the mass spectrum at generator level is 
not distorted by this reweighting procedure. Figure 10 shows the linear dependence of Amw 
on the average kaon and proton multiplicity. This linear dependence is used to translate the 
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uncertainty of the measured kaon and proton multiplicities into uncertainties on mw and Fw 
Table 10 presents the shifts and uncertainties of mw and Fw due to the correction of the Monte 
Carlo simulation. 

In the qqqq final state, the DURHAM parameter 7/45 is used to discriminate events with 
hard-gluon radiation. This variable might be affected by hadronisation uncertainties. A change 
of the selection criterion logy45 < —6.2 between —5.8 and —6.6 has no significant infiuence on 
mw 

The average absolute shift of m-w and Fw due to the alternative hadronisation models 
ARIADNE and HERWIG and the uncertainty deduced from the variation of the kaon and 
proton multiplicities in the Monte Carlo simulation are added in quadrature to yield the total 
uncertainties due to the hadronisation modelling, given in Tables 4 and 5. 

6.9 Final state interactions in fuUy-hadronic events 

The Monte Carlo programs hadronise the quarks from the two W bosons independently. How- 
ever, CR effects would invalidate this assumption and thus affect the mass reconstruction. 
Similarly, BEC between bosons arising from different W bosons, if incorrectly modelled, could 

have the same effect. 

Our measurements [22,39,40] of BEC indicate that correlations in hadronic W-boson decays 
are very similar to those in Z-boson decays into light quarks. Furthermore, BEC between 
hadrons from different W bosons are disfavoured. They are limited to at most 30% of the 
strength simulated in the BE32 model [21] implemented in PYTHIA 5.7. Since all our previous 
mass measurements at ^/s = 172 183 GcV were performed under the assumption of full inter- 
W BEC, the results obtained in the qqqq channel are re-evaluated in light of our measurement 
of vanishing inter- W BEC. 

Reference 22 presents the L3 measurement of the difference between the two-particle den- 
sities of the data and the simulation without inter-W BEC, Ap2{Q). The integral, J, of this 
difference is measured to be below 0.39 at 68% confidence level. For different Monte Carlo 
samples, generated with various strengths of inter-W BEC, but fixed strength of the intra-W 
BEC, the integral J is determined. The shift Amw exhibits a linear dependence with respect 
to J, as shown in Figure 11. The effects on mw and Fw for a maximum inter-W BEC, as 
allowed by our direct BEC measurement, are detailed in Table 11. A linear dependence of BEC 
effects on ^/s is assumed. 

A dedicated study of reconnection effects in the particle fiow between jets in qqqq events 
shows that the data are consistent with no or only a small CR effect [24] . A 68% upper limit on 
the CR parameter kj is set at 1.1 in the framework of the SK-I model [41] as implemented in 
PYTHIA 5.7. The influence of the CR parameter ki on mw is studied in the SK-I framework by 
mixing event samples simulated at ^/s = 189 GcV with full and without inter-W CR. The result 
is shown in Figure 12 for moderate values of ki where a linear dependence can be assumed. 
The particle flow analysis is found to be insensitive to CR effects implemented in other models 
such as ARIADNE type II [42] and HERWIG [43]. The ARIADNE-II model is compared to 
the ARIADNE-I model, the latter having been modified such that in both models the shower 
cascade is performed in two phases with an identical cut-off parameter. 

These Monte Carlo studies show that the effect of CR on mw grows with increasing ^/s 
in the case of the SK-I model, while only httle dependence on ^/s is seen for ARIADNE and 
HERWIG. For all energies and all models the shift of mw is comparable or smaller than the shift 
predicted by the SK-I model at ki — 1.1, which is used to estimate the systematic uncertainty 
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due to CR effects. It is interesting to note that studies of the distribution of particles in the 
inter-jet region of three-jet hadronic Z decays exchidc the predictions of the CR models of 
ARIADNE and HERWIG for this case [44]. No version of the SK-I model apphcable to Z 
decays exists. 

The use of a cone algorithm for jet clustering lowers the sensitivity to CR effects, as the 
analysis is less affected by the inter-jet regions where the influence of CR is largest. More effec- 
tively, removing clusters below a certain energy cut rejects particles predominantly produced 
during the non-perturbative phase of the hadronisation process where CR effects take place. 
Monte Carlo studies are performed at \/i = 189 GeV applying various cuts on the minimum 
cluster energy. The dependence of the mw shift on the energy cut is extrapolated to the full 
data sample and shown in Figure 13. The additional component to the statistical uncertainty 
due to the slight degradation of the mass resolution caused by the cut is calculated and added 
in quadrature to the shift of mw- A cut at minimum cluster energies of 2 GeV is found to be 
the optimal choice and is therefore used in the extraction of mw and Fw from the data of the 
qqqq final state. Table 12 presents the effect of CR on mw and Fw- 

Monte Carlo studies show that the relative reduction of the mw shift due to the energy cut 
is independent of ki and ^/s. The mass shifts observed for the SK-I Monte Carlo simulation 
with full CR at various ^/s values are obtained using the dependence on ki and on the energy 
cut extracted at ^/s = 189 GeV. The systematic uncertainties on mw are calculated using a 
linear dependence on ^/s and assumed to be fully correlated. For Fw no ^/s dependence is 
seen. 



Figure 14 compares the mw measurements in the four different final states at the seven average 
■\/s values. The measurements of mw and Fw from the individual final states are combined 
using the "best hnear unbiased estimate" technique [45]. This combination method takes into 
account all systematic uncertainties as well as their correlations. When combining measure- 
ments taken at different ^/s values, the correlations of the LEP energy determination [27] are 
used. Within each final state, the uncertainties due to Icpton measurement, background de- 
termination, BEG and CR are taken as fully correlated between the measurements at different 
^/s. The uncertainties due to jet measurement, photon radiation and hadronisation are fully 
correlated between all final states and between all ^/s values. The systematic uncertainty due 
to limited Monte Carlo statistics remains uncorrelated for all measurements. In the case of the 
simultaneous estimate of mw and Fw, the correlations between both parameters as determined 
in the individual box fits are included in the combination procedure. 

Combined results of mw are shown in Figure 15 for each y/s value averaged over the final 
states. Figure 16 shows the results for each final state and their combination. Table 13 gives 
the results on mw for each final state. The combination of the results at ^/s — 189 — 209 GeV 
yields for the semi-leptonic and the fuUy-hadronic final states: 



Here and in the following, the first uncertainty is statistical and the second systematic. The 
qqlv and the qqqq channels exhibit a correlation of 9%. The contributions of the individual 
sources of systematic uncertainty to the combined mw value in the qqiu channel is given in 
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mw(qq^i^) = 80.196 ± 0.070 ± 0.026 GeV and 
^w(qqqq) = 80.298 ± 0.064 ± 0.049 GeV . 



(8) 
(9) 



Table 4. 
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The difference between the values of mw determined in the qq£u and qqqq channels is 

mw(qq^i^) - mw(qqqq) = -0.088 ± 0.094 ± 0.031 GeV . (10) 

BEG and CR effects are not included in the systematic uncertainty on the mass difference. 
Moreover, hadronisation uncertainties are treated as uncorrelated between the qqiu and qqqq 
final states. This causes the mass difference not to equal the difference of the mass values given 
in Equations (8) and (9). 

Averaging the values of the qqiu and qqqq channels, including BEG and GR uncertainties 
and all correlations, yields 

mwiffff) = 80.242 ± 0.048 ± 0.031 GeV . (11) 

In this combination the value of x^/d.o.f. is 29.2/27 and the weight of the fuUy-hadronic 
channel is 46%. In absence of any systematic uncertainties, the statistical precision of the 
measurement would be 47 MeV. In Table 4 the contributions of the individual sources of 
systematic uncertainty to this combined mw result are given. 

The results in this paper are combined with the direct measurements obtained at y/s — 
172 - 184 GeV [7] to give 

mw(qqM = 80.212 ± 0.066 ± 0.027 GeV and (12) 
mw(qqqq) = 80.325 ± 0.061 ± 0.052 GeV , (13) 

with a correlation of 10%. Combining the results from direct measurements at ^/s = 172 — 
209 GeV with those result obtained from cross section measurements at ^/s = 161 — 172 GeV [5] 
yields 

mw = 80.270 ± 0.046 ± 0.031 GeV . (14) 

The W-boson width is determined in fits for both mw and Fw Table 14 gives the results 
for ^ = 189 - 209 GeV. Gombining all data yields 

Tw = 2.18 ±0.11 ±0.09 GeV . (15) 
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Table 1: Integrated luminosity, £, together with the number of selected data events, iVdata, and 
expected number of background events, A^back, for each final state and average value of ^/s. 





Eiiorg}- ['X] 


[clog.] 


o [dog.] 


Electrons 


1.4 


0.47 


0.083 


Muons 


5.2 


0.22 


0.007 


Hadronic jets (no cut) 


15 


2.4 


1.9 


Ha(h-()iiic jets {Ec > 2 GeV) 


15 


2.5 


2.1 



Table 2: Average energy and angle resolutions for reconstructed electrons, muons and hadronic 
jets as determined in Monte Carlo simulation. Resolutions for hadronic jets are given with and 
without the cut on the minimum cluster energy, Ec- 
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Table 3: Mass resolutions in GeV as determined in Monte Carlo simulation: raw mass resolution, 
m^'aw) of the hadronically-decaying W bosons; resolution of the average of the two raw masses in 
each event, m^w; resolution after rescaling the jet energies, mresc, or after applying a kinematic 
fit, m„c- The last column indicates which of the mass variables is used for each final state in 
the extraction of mw and Fw, as described in Section 5. 
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Total statistical 


99 


119 


175 


64 


70 


48 


Total 
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74 


57 



Table 4: Systematic uncertainties on mw, in MeV, for the various final states. The values refer 
to the complete data set at y/s = 189 — 209 GeV and take into account correlations between 
energy points and final states. 
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65 


90 


180 


115 


85 


90 


Total statistical 
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Total 


255 


315 


420 
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Table 5: Systematic uncertainties on in MeV, for the various final states. All uncertainties 
are rounded to the next 5 MeV. The values refer to the complete data set at ^/s — 189—209 GeV 
and take into account correlations between energy points and final states. 



21 







Amw 


[MeV] 




AFwl 


[MeV] 


qqeu 


qqjiu 


qqriy 


qqqq 


qqeu 




qqriy 


qqqq 


Electron energy 


6 










12 










Muon energy 






12 










37 






Jet energy scale (±50 McV) 


3 




10 


9 


2 


1 




9 


16 


4 


Jet energy smearing (1%) 


1 




4 


7 


1 


10 




25 


53 


7 


Jet angle smearing (0.5°) 


2 




4 


20 


4 


17 




16 


47 


18 



Table 6: Changes of mw and Fw due to variations of the energy measurement of electrons, 
muons and jets and the resolutions of the jet directions. 
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Table 7: Changes of mw and Fw due to variations of the background processes. For fuUy- 
hadronic events the uncertainty due to the yu spectrum is also given. 
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Monte Carlo reweighting 
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Table 8: Changes of mw and Fw due to variations in the modelling of photon radiation. The 
first row gives the results of a comparison between the RACOONWW and KANDY generators. 
The second row gives the difference between the 0{a^) and the 0{a'^) calculation, obtained 
by reweighting KANDY events. The statistical accuracy of the generator comparison is about 
8 MeV, while the statistical uncertainty of the reweighting procedure is negligible. 
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Amw [MeV] 


AFw 


MeV] 


qqev 


qqjiu 


qqri/ 


qqqq 


qqeu 


qqjiiy 


qqru 


qqqq 


HERWIG 
ARIADNE 




-15 


-8 
-11 


-41 
-44 


-3 

+11 


-96 
+15 


-141 
-1 


-275 
-24 


-32 

+5 



Table 9: Changes of mw and Fw due to the use of the hadronisation models HERWIG and 
ARIADNE instead of PYTHIA. The statistical accuracy is always better than 15 MeV and 
30 MeV for mw and Fw, respectively. 





Amw [MeV] 


AFw [MeV] 


qqiu 


qqqq 


qqiu 


qqqq 


Kaon multiplicity 
Proton multiplicity 


+13 + 12 
+1+2 


+25 ± 23 
+7+15 


-12 ± 11 
-3+5 


-95 + 87 
-36 + 80 



Table 10: Changes of mw and Fw due to reweighting Monte Carlo events with respect to 
variations of the mean charged-kaon and proton multiplicities. The given uncertainties are due 
to the experimental uncertainties in the determination of these multiplicities [38] . 



[GeV] 


Amw [MeV] 


AFw [MeV] 


189 
207 


+11 
+23 


-1 
-15 



Table 11: Changes of mw and Fw in the qqqq channel when replacing our standard simulation 
by the PYTHIA BE32 model with a strength of inter-W BEC corresponding to the 68% upper 
limit set by our direct BEC measurement [22]. The statistical accuracy is 6 MeV and 14 MeV 
for mw and Fw, respectively. 
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Amw [MeV] 


AFw [MeV] 


189 
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-22 

-57 


-48 
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Table 12: Changes of mw and Fw in the qqqq final state when replacing our standard simulation 
by the PYTHIA SK-I model with h = 1.1 which is the 68% upper limit set by our CR 
measurement [24] . The cut on the minimum cluster energy of 2 GeV is apphed. The statistical 
accuracies are about 10 MeV for mw and 20 MeV for Fw- 
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Process 


mw [GeV] 


<S [GeV] 


e''"e~ — > qqeu 
e'^e~ — > qqptiy 
e'^e~ — > qqriy 


80.225 ± 0.099 ± 0.024 
80.152 ±0.119 ±0.026 
80.195 ±0.175 ±0.060 


0.095 
0.119 
0.162 


e'^e~ — > qqiu 


80.196 ±0.070 ±0.026 


0.068 


e'^e~ — > qqqq 


80.298 ± 0.064 ± 0.049 


0.062 


e+e- ^ //// 


80.242 ± 0.048 ± 0.031 


0.047 



Table 13: Results on mw for the data collected at ^/s = 189 — 209 GeV. The first uncertainty is 
statistical and the second systematic. Also shown is the expected statistical uncertainty, (7^^^^. 



Process 


mw [GeV] 


Fw [GcV] 


Correlation 


e"'"e~ qqiu 
e"'"e~ — qqqq 


80.174 ±0.078 ±0.027 
80.284 ±0.074 ±0.050 


2.50 ±0.17 ±0.09 
1.97± 0.15 ±0.12 


0.01 
0.15 


e+e- ^ //// 


80.236 ± 0.054 ± 0.032 


2.22 ±0.11 ±0.09 


0.14 



Table 14: Results on mw and Fw obtained from a simultaneous fit of both quantities using 
data collected at = 189 — 209 GeV. The first uncertainty is statistical and the second 
systematic. Also quoted is the correlation between mw and Fw- 
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Figure 1: Improvement of mass resolutions due to kinematic constraints for a) qqez/, b) qq/ii^, 
c) qqru and d) qqqq events. The open circles represent the raw mass spectra and the full points 
the spectra obtained after applying the kinematic fit or the jet-energy rescaling. Monte Carlo 
predictions are also shown. In a), b) and d) m^aw is the average of the two raw masses while in 
c) m'^2w mass of the hadronic system. 
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Figure 2: Distributions of reconstructed W-boson masses after applying the kinematic fit using 
the equal-mass constraint for the a) qqez/, b) qq/iz/ and c) qqru channels and d) the best pairing 
for the qqqq channel. The signal Monte Carlo events are reweighted according to the fitted 
value of mw- 
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Figure 3: Distribution of reconstructed W-boson masses after applying the kinematic fit using 
the equal-mass constraint for semi-leptonic final states. The signal Monte Carlo events are 
reweighted according to the fitted value of mw- 
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Figure 4: Distribution of reconstructed W-boson masses after applying the kinematic fit using 
the equal-mass constraint for all W pairs. The signal Monte Carlo events are reweighted 
according to the fitted value of mw- 
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Figure 5: Calorimetric energy-flow versus the angle relative to the direction of the charged 
lepton, (, for a) the qqeu and b) the qq/iu events. The error bars represent the standard 
deviation of the data distribution in each bin. 
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Figure 6: Distribution of the mass shifts between the standard analysis of qqqq events and 
analyses using a) a displaced vertex and b) jet reconstruction from tracking information only. 
A Gaussian fit is applied to the data distribution of a) and indicates an average mass shift 
consistent with zero, as shown by the curve. The data distribution of b) is in good agreement 
with the Monte Carlo prediction. 
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Figure 7: a) Energy spectrum of the clusters used in the jet reconstruction and changes of mw 
for b) the qqiu and c) the qqqq final states caused by a variation of the cut on the minimum 
cluster energy, Ec- The arrows show the default values of the cut. 
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Figure 8: a) Energy flow as a function of the angle relative to the jet direction, and changes 
of mw for b) the qqiu and c) the qqqq final states after removing clusters outside a cone of 
half-opening angle ( around the jet direction. 
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Figure 9: a) Comparison of the reconstructed mass spectra, combined for all final states at 
a/s = 189 GeV, for data and for the three hadronisation models PYTHIA, ARIADNE and 
HERWIG and b) the spectra normalised to the PYTHIA expectation. 
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Figure 10: Changes of mw due to reweighting Monte Carlo events according to the mean 
charged-kaon multiphcity for a) the qqii^ and b) the qqqq events and of the mean proton 
multiphcity for c) the qq£u and d) the qqqq events. The full circles show the default values of 
our simulation whereas the vertical lines show the measured multiplicities and the grey bands 
their uncertainties [38]. 
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Figure 11: Changes of m-yv with respect to the observable J for Monte Carlo samples of the 
BE32 model with different BEC parameters at ^/s = 189 GeV. The grey band shows the range 
of J which is compatible with our BEC measurement [22] at the 68% confidence level. 




Figure 12: Changes of mw with respect to the parameter kj of the SK-I model at ^/s = 189 GeV. 
The cut on the minimum cluster energy of 2 GeV is applied. The grey band shows the range 
of ki which is compatible with our CR measurement [24] at the 68% confidence level. 
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Figure 13: CR effects simulated witli tlie Monte Carlo model SK-I calculated after removing 
clusters with an energy below a given threshold energy. The change of the final m-yv measure- 
ment in the qqqq channel, Amgyst, when the default simulation without CR is replaced by the 
SK-I model using ki = 1.1 is shown. The additional component of the statistical uncertainty 
on the final mw result, Acistat, after applying the given energy cut and the quadratic sum of 
both effects is also shown. 
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Figure 14: The results for m-yv for the four final states and the seven average ^/s values. 
Statistical and systematic uncertainties are added in quadrature. The combined mw result and 
its uncertainty are indicated as the dashed line and the grey band, respectively. 
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Figure 15: Comparison of the results for for the seven average ^/s values. The inner error 
bar represents the statistical uncertainty. The combined mw result and its uncertainty are 
indicated as the dashed line and the grey band, respectively. 
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Figure 16: Comparison of the results for mw for the = 189 — 209 GeV in each of the different 
final states. The inner error bar represents the statistical uncertainty. 
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